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Abstract. I describe how the gluon chain model of QCD string formation meets a 
number of criteria which are required of any theory of the confining force, includ- 
ing: the correct center dependence and (at large-N) Casimir scaling of the string 
tension, the logarithmic broadening of the QCD flux tube, and the existence of a 
Liischer term in the static quark potential. 



1. Introduction 

The gluon chain model [1-4] is a picture of the formation and composition of the 
QCD string in terms of the perturbative excitations of the theory. In this talk I 
would like to discuss some recent work [5] in this area, which was carried out in 
collaboration with Charles Thorn. 

Normally gluon (particle) excitations are useful for describing high-energy 
scattering processes, while non-perturbative effects, such as confinement and chi- 
ral symmetry breaking, are usually ascribed to some special class of field config- 
urations with particular topological properties. However, it is not excluded that 
particle and field descriptions of various phenomena in QCD can overlap in some 
cases. An example is the color screening of higher representation Wilson loops. 
In this case one is easily convinced, just by thinking about string breaking due to 
gluon pair production, that the asymptotic string tension of a higher-representation 
loop can depend only on the transformation properties of the representation with 
respect to the center subgroup (i.e. on the "N-ality" of the representation). For an 
adjoint string in particular, one ends up with two gluelump states, each consisting 
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of a gluon bound to a heavy source. On the other hand, from the "field" point of 
view, the area law falloff of the Wilson loop is due to large scale vacuum fluctua- 
tions of some kind. The N-ality dependence of the string tension tells us that these 
large-scale field fluctuations must induce fluctuations in Wilson loop holonomies 
only among the center elements, rather than coset elements, of the SU(N) gauge 
group. This simple fact implies (very strongly, in my opinion) the existence of a 
center vortex mechanism of some kind. 

In principle we then have two complementary ways of understanding the N- 
ality dependence of the QCD string tension: (i) in terms of color-screening by 
constituent gluons; and (ii) in terms of the vortex confinement mechanism. But if 
we can picture string breaking in terms of particle excitations, is it also possible to 
describe string formation in terms of constituent gluons? As it happens, large-A^. 
considerations lead to just such a description. Consider a very high-order planar 
Feynman diagram contributing to a Wilson loop expectation value, as shown in 
Fig. 1. A time-slice of the diagram reveals a sequence of gluons, interacting only 
with their nearest neighbors in the diagram. If we take this picture seriously, it 
suggests that the QCD string might be regarded, in some gauge, as a "chain" of 
constituent gluons, with each gluon held in place by its attraction to its two nearest 
neighbors in the chain. 



timeslice 
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Figure 1. The gluon chain as a time slice of a planar diagram (shown here in double-line notation). 
A solid (open) hemisphere indicates a quark (antiquark) color index. 



The linear potential in this "gluon chain" model comes about in the following 
way: As the heavy quarks separate, we expect that at some point the interaction 
energy increases rapidly due to the running coupling. Eventually, it becomes ener- 
getically favorable to insert a gluon between the quarks, to reduce the color charge 
separation. As the quarks continue to separate, the process repeats, as shown in 
Fig. 2, and we end up with a chain of gluons. The average gluon separation R 
along the axis joining the quarks is fixed, regardless of the quark separation L, 
and the total energy of the chain is the energy per gluon times the number N of 
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gluons in the chain, i.e. 

Echain ~ NE g i uon = g ^°" L = CL (1) 

where E g i uon is the (kinetic+interaction) energy per gluon, and a = E g i uon /R is the 
string tension. In this picture, the linear growth in the number of gluons in the 
chain is the origin of the linear potential. 

in. - • * -o » o- 
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Figure 2. Constituent gluons (small circles) appear as quarks separate, to keep the average color 
charge separation below some maximum value. Dotted lines indicate nearest-neighbor interactions. 

Every theory of confinement aims at explaining the linearity of the static quark 
potential at large distances. However, it is by now well understood that linearity 
is only one of several conditions that a theory of the confining force must satisfy. 
The requirements also include: 

1. Casimir Scaling. At intermediate distance scales, the string tension should 
be proportional to the quadratic Casimir of the quark color group representa- 
tion. 

2. Center Dependence. Asymptotically, the string tension should depend only 
on the N-ality, not the Casimir, of the quark color group representation. 

3. String Behavior: 

— Roughening. There is a logarithmic broadening of the flux tube with 
quark separation L. 

— Liischer term. There is a universal —c/L term in the static quark po- 
tential at large distances. 

Taken together, this is a challenging set of conditions. The abelian projection 
theory, for example, has difficulties with both Casimir scaling [10] and center de- 
pendence [11]. Other proposals, which seek to derive confinement from a. \/k 4 
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behavior of the gluon propagator in some gauge, will run into trouble when con- 
fronted with roughening and the Liischer term. So does the gluon chain model, 
which is yet another idea about confinement, do any better in meeting these con- 
ditions? 

Having already discussed the linearity of the static potential, let us consider 
Casimir scaling. In any S\J(N C ) representation r, the group character can be ex- 
pressed as a product 

lr[g] « {%F[g]) (x* F [g}) + sub-leading terms. (2) 

where F denotes the fundamental representation. By factorization at large-A^., a 
Wilson loop in representation r has a string tension 

c r = M r a F (3) 

at N c — > oo, where M r = n + n. In this limit, the quadratic Casimir is C r = M r N c /2. 
Exact Casimir scaling is therefore a property of the planar limit. The gluon chain 
model, which is motivated by large-A^. considerations, inherits this property. In 
the gluon-chain model, there are M r chains terminating at each heavy source, 
and these are non-interacting at N c —> °°. The total energy of the system is then 
proportional to the number of chains times the length of each chain, and in this 
way Casimir scaling 

O r oc C r {N c -► oo) (4) 

is obtained at large N c . The situation for heavy adjoint-representation quarks, 
which have two chains between them, is shown in Fig. 3(1). 

The appropriate center dependence of the string tension is due to color screen- 
ing. Screening is accomplished by a (1/Af? suppressed) contact interaction be- 
tween constituent gluons in different chains, leading to string-breaking processes 
such as those shown in Fig. 3. 

As for string behavior, the gluon chain is clearly some sort of discretized 
string, so string-like properties are not entirely unexpected. I will return to this 
topic in section 4, below. 



2. The Force Renormalization Scheme 

The growth of the running coupling with color charge separation is an essential 
ingredient of gluon chain picture, but the relationship of the running coupling to 
the interaction potential, which depends on the choice of renormalization scheme, 
is an important consideration. Suppose one defines the running coupling in terms 
of the static potential via 
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Figure 3. Adjoint string-breaking in the gluon chain model. Two gluons in separate chains (I) 
scatter by a contact interaction, resulting in the re-arrangement of color indices indicated in II. This 
corresponds to chains starting and ending on the same heavy source. The chains then contract down 
to smaller "gluelumps" (III). 



Actual 




running coupling? 



Figure 4. A running coupling in the V-scheme which is monotonically increasing takes the potential 
in the negative-V direction, away from the actual potential. 



This is the V-renormalization scheme. But if grows monotonically, the result- 
ing behavior of V (R) is the opposite of what is required for confinement, behaving 
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roughly like the curve labeled "running coupling" in the sketch shown in Fig. (4). 
In order that the potential becomes positive at some point, as it does in the actual 
potential found in numerical simulations, it is obviously necessary that a% (R) be- 
comes smaller and eventually changes sign as R increases. This is asking a lot of 
perturbation theory, which tends in the opposite direction. A better renormaliza- 
tion scheme, advocated by Grunberg [7], Sommer [8], and also by Thorn, is to 
define the running coupling in terms of the force 

and derive the potential at intermediate distances from integration 



V 9 g(R)=V 9 g(R A )+ j R dR 

JRa 



F(R) 



(7) 



In this "F-scheme" the running coupling doesn't have to change sign, and in fact 
the three loop term in the beta function 

3 / 11 102 2 1.65 4 \ 

is substantially smaller (by a factor of 2.6) in the F-scheme than the V-scheme. The 
three-loop potential has been computed in the F-scheme by Necco and Sommer 
in ref. [6]. I have extracted Fig. 5 below from their article; this figure displays 
the perturbative potential compared to the corresponding Monte Carlo data. The 
parameter r<) ~ 0.5 fm is the Sommer scale. Notice that the F-scheme gives a 
remarkably good fit to the numerical data almost up to the Landau pole, while the 
V-scheme results in a very poor fit to the data in this interval. 

According to a rigorous theorem [9], the static potential must be concave 
downwards; i.e. the force cannot increase with distance. This implies that the 
the perturbative 3-loop result must certainly break down at quark separation L m 
0.6ro « 0.3 fm in the F-scheme. In the gluon-chain model, the rising force is 
avoided by inserting constituent gluons between the quarks. From the breakdown 
of perturbation theory at L = 0.6ro, we can estimate that the inter-gluon separation 
along the line joining the quarks is R = 0.3ro. 

3. A Variational Approach 

Let *Fo [A] be the QCD vacuum wavefunctional. An excited state can be written as 

V„[A]=fi[A]V [A]. (9) 

e.g. for a glueball 

*g[A]= £vg°EA] (10) 

N=\ 
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Figure 5. The static potential in the F and V-renormalization schemes, compared to numerical data. 
This figure is taken from Necco and Sommer, ref. [6]. 



where 



Wq^IA] = dx\dx2-..dxN Ui ltl ...i lN {xuX2,...,x N ) 

TrA w (x { )A K (x 2 )...A m (x N ) \ ¥ [A] 



E = 



is the N constituent gluon state. The energy is 

Q t = Q[A(x,t) 



(H) 



Defining 
it is easy to show that 



-- lim — log(oto_ r ) 



(11) 

(12) 
(13) 

(14) 
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where (...) represents the Euclidean VEV. 

The idea is to choose Q as a function of a few parameters, and then determine 
the optimal parameters by minimizing the value of E computed perturbatively [3]. 
For the gluon chain, we propose to use 

Qchain [A] = q ai (0){J dx\ dx 2 . . .dx N y,t m ... m (xi,X2,...,X N ) 

A^(xOA£f (x 2 )..A™ (x N )}r N+> (L) (15) 

Before trying this approach in the full-blown field theory, it is interesting to 
apply these ideas to a simple quantum-mechanical model with most of the same 
qualitative features. Let us take the multi-gluon Hamiltonian to be 

N-l N-l 

H= £ P„ + £ ViXn-Xn-J+Vqgixti+VqgtL-XN-!) (16) 

n=l n=2 

The Product Ansatz for the (N-l)-gluon chain ground state is 

N 

*(xi,X2,...,^ v _i)=Af[\|/(u ! ) (17) 

1=1 

where 

u i =x i -x i -i (18) 

and 

¥ (u) = e - ulllr \ (19) 

with the constraint 

xq = , xn = L (20) 
After some work, one finds that the string tension (= energy/unit length) is 

!=^+>»» <21 > 

In particular, taking for V the instantaneous Coulomb potential V(u) = — CfOL s /\u\, 
we have 

£ 1 ^ C ^erf(^l. (22) 



L rR\K R 2 



where both r and N (no. of gluons), or equivalently r and R = L/N, are taken as 
variational parameters. 

Our second trial wavefunction, the String Ansatz, is adapted from discretized 
light-cone string theory [12]: 

a i m *i>(x 1 ,x 2 ,...,x N - 1 )=0 | ^2y N ~ l (23) 
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with 



x k 



Pk 




1 t-. -t\ • 



T L \J^-{ a m-<)^{ -rrk) (24) 

m=\ ' 



where co m = 2 sin S^. Define r = 2/ ^/STo- For a Coulomb potential, we then find 



1 8 C F a 



L rR 7I 3 /2 R2 



(25) 



Comparison with the product ansatz shows that the string ansatz achieves a slightly 
lower energy. 

To estimate the string tension, we begin from 



£ _ 1 
T~R 



5^7 + ("CO 



and make the approximation that (V(u)) = V(s), where 

~s = (w\) = R eri (j : ) 

Minimizing wrt R and r results in two conditions 



(26) 



(27) 



£ 
L 
1 



n 3/2 r 2 



which together imply 



ds_ 
dR 
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(28) 



(29) 



Given F(s) = dV /ds, we could determine r,R and a. Even without F(s), we 
can still make a "ballpark" estimate of a. First, the static force increases beyond 
L = 0.6ro, so a constituent gluon must appear between the quarks at that point, 
to stabilize the force. This suggests that R = 0.3ro. Secondly, the running cou- 
pling a(s) must be 0(1), if the growth in potential energy over-compensates the 
falloff in kinetic. This argues for s « 0.6ro and r « 0.63ro. With these guesses, 
one finds a = 1.26/ r^, which is to be compared with the experimental value 
Gexp = (430 Mev) 2 = 1.18/7-q. We conclude that the value of the string tension 
may work out in a genuine field-theory/variational calculation, but of course this 
remains to be seen. 
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4. String Behavior 



Using the string ansatz wavefunction, one rinds for the transverse coordinate of 
the k-th gluon in the chain 



<o|4i|o> 



D-2 



iV-l 



sin 2 (milk/ N) 
IMh^ sin(mn/2N) 



■InN 



D-2. 
2nT 

r 2 {D-2) 



In 



R 



(30) 



which demonstrates the expected logarithmic broadening of the flux tube with 
quark separation. Note that there is no need to insert a high-frequency cutoff; this 
is taken care of by the discreteness of the gluon chain. 

The total energy of the gluon chain in the string ansatz can also be computed, 
and the result minimized wrt the variational parameters. The result is 

■E = ^ 2 Jd\e-» 2 l.Vv(R + 2u^) 



(31) 







( 4R )- 




\rs/Z) 





rrc 3 / 2 



(see ref. [5] for details). The first term on the rhs is the linear potential, the second 
is a Liischer-like term. For the bosonic string, the factor multiplying n/24L in 
the Liischer term would be D — 2 = 2. With our estimate above of R r/2 this 
coefficient is estimated to be 2/y / 7TRi 1.13. At this stage, however, numerical 
values should not be taken too seriously. 

5. Numerical Studies 

It is possible to investigate the constituent gluon composition of the QCD string 
by numerical methods, and in fact a study along these lines was carried out many 
years ago (second reference of ref. [2]). The idea is the following: Observe that 
the ground state of the QCD flux tube, in a physical gauge, can be written as 

^stringlAfifat = 0)] = 

lim (w(R,2T)Z^j 1/2 J DA^{x,t < 0)8[F(A)]A[A] 

xPexp[ij) rfjc%Jexp - J dt L[A] (32) 



11 



t = 



t = -T 



Figure 6. Contour for generating the QCD string state. 



where C is the R x T open contour shown in Fig. 6. The R x 2T Wilson loop 
factor W(R,2T) ensures that ^string is normalized to unity. 

If the gluon chain model is correct, then ^string can a l so be expressed as a sum 
of n-gluon states 



string 



¥„[A] 

where Q n is an n-gluon operator, and 

,1 



£c n ¥„[A] 



N n = {-Tr[QlQ n ])-" 2 



(33) 



(34) 



is a normalization constant. Assuming the *P„ states are orthogonal, the overlap of 
an n-gluon state with the true string state is given, in the T — > °° limit, by 



Tr[2jPexp[/<f c _dx?A f . 



Tr[QlQ n ])W(R,2T 



(35) 



The second article of ref. [2] used the following trial operators on the lattice 
(with A„ = (U fl -Ul)/2i); 



Qo 



ab 



Qf 



tab 



I Af (x) 

X\ <X<X2 

Xl<X<X2X<X t <X2 



where 



A a x h (x)=£A?(x,y,z)e 



-5r, 



(36) 



(37) 
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is a "smeared" vector potential, rj_ is the transverse distance from the flux tube 
axis, and 8 is a variational parameter. The corresponding zero, one, and two-gluon 
overlaps (cq, c\, ci) were then computed from eq. (35) via lattice Monte Carlo. 
It was found that for small R, the zero-gluon overlap cq is the largest of the three 
overlaps. As R increases, the Co overlap falls and the c\ term becomes dominant, 
and at the largest R separations that were used in the simulations we found the two- 
gluon overlap C2 becoming the largest of the three. The sum of the three overlaps 
accounts for about 90% of the norm of ^string, which means that for the given 
range of R the sum of zero, one and two gluon states is a fairly good approximation 
to the true flux tube state. For details, please see the cited reference. 

All of this is in good qualitative agreement with the gluon chain picture, but 
the old results can surely be much improved. It may be useful to repeat the inves- 
tigation of ref. [2] with modern computers and better noise-reduction techniques. 

6. Conclusions 

The gluon-chain model offers a simple and concise explanation of many features 
of the confining force — Casimir scaling at large-/V c , center dependence, rough- 
ening, and the Liischer term — which are problematic in many other approachs. 
The model is essentially a "particle" picture of string formation, and I regard it as 
complementary to the "field" explanation of confinement in terms of center vor- 
tices. The next step will be to apply the field theory/variational approach outlined 
above to obtain quantitative estimates for the string tension and, perhaps, for the 
masses of the low-lying glueballs. 

References 

1. C. B. Thorn, Phys. Rev. D19 (1979) 639; D20 (1979) 1435; D20 (1979) 1934. 

2. J. Greensite, Nucl. Phys. B249 (1985) 263; B315 (1989) 663. 

3. J. Greensite and M. Halpern, Nucl. Phys. B271 (1986) 379. 

4. C. B. Thorn, Phys. Rev. D51 (1995) 647. 

5. J. Greensite and C. B. Thorn, J. High Energy Phys. 02 (2002) 014; hep-ph/01 12326. 

6. S. Necco and R. Sommer, Phys. Lett. B523 (2001) 135; hep-ph/01 09093. 

7. G. Grunberg, Phys. Rev. D40 (1989) 680. 

8. R. Sommer, Nucl. Phys. B411 (1994) 839. 

9. C. Bachas, Phys. Rev. D33 (1986) 2723. 

10. L. Del Debbio, M. Faber, J. Greensite, and S. Olejnik, Nucl. Phys. Proc. Suppl. 53 (1997) 141, 
hep-lat/9607053. 

11. J. Ambj0rn, J. Giedt, and J. Greensite, J. High Energy Phys. 02 (2000) 033, hep-lat/9907021. 

12. C.B. Thorn, Nucl. Phys. B263 (1986) 493. 



